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I.  INTRODUCTION 


Spike-nosed  projectile  configurations  are  used  against  armored  targets  where  the  spike 
provides  a  stand-ofF  distance  between  the  armor  and  the  shaped-charge  warhead.  These 
projectiles  have  been  continuously  developed  and  improved,  despite  their  high  drag  disad¬ 
vantage.  One  major  difficulty  with  the  spike-nosed  configurations  is  the  possibility  of  two 
flow  patterns.  In  many  cases,  one  does  not  know  which  flow  pattern  will  occur  under  what 
conditions.  One  flew  pattern  is  of  low  drag  and  the  other  is  of  a  higher  drag.  Unfortunately, 
the  high  drag  pattern  is  more  dominant  and  more  stable  while  the  lower  drag  pattern  is  less 
dominant  and  occurs  in  fewer  cases.  Some  projectiles  launched  at  high  supersonic  speeds 
start  with  the  high  drag  and  then  switch  to  the  lower  drag  mode  which  is  more  stable  at  a 
lower  speed.  In  other  cases  the  flow  may  oscillate  slowly  between  the  two  modes.  In  other 
configurations  the  flow  oscillates  rapidly  causing  the  “buzzing”  phenomenon.  The  critical 
Mach  number  for  switching  patterns  depends  on  the  Mach  number,  length  of  the  spike,  spike 
diameter  to  projectile  diameter  ratio,  and  to  a  lesser  degree  on  the  Reynolds  number  and 
flow  turbulence. 

Platou*  reported  the  two  flow  patterns  experienced  during  supersonic  wind  tunnel  test¬ 
ing  in  the  early  1950’s.  Koenig  et  aU  more  recently  reported  hysteresis  mode  changes  during 
transonic  wind  tuimel  testing.  Shang  and  Hankey^  successfully  and  accurately  simulated 
the  buzziiig  phenomenon  in  the  1970’s  for  a  re-entry  body  with  a  relatively  short  spike. 
Calarese  and  Hankey'*  tested  and  analyzed  the  oscillations  of  spike-tipped  bodies.  Mikhail® 
in  1989  successfully  computed  the  flow  past  three  sharp-edged,  spike-nosed  configurations 
and  induced  the  two  modes  for  each  case.  Mikhail  verified  his  results  at  Mach  =  1.72  with 
those  wind  tunnel  tests  of  Reference  1. 

To  avoid  the  uncertainity  of  the  dual  flow  modes,  it  had  been  found  long  ago  that 
by  inserting  a  ring  near  the  spike  tip  the  high  drag  mode  is  prevented  by  introducing  flow 
separation  along  the  spike  and  preventing  reattachment  before  the  flow  faces  the  projectile 
shoulder  which  causes  a  shoulder  “bow”  shock  and  the  high  drag.  No  specific  name  is 
credited  in  the  literature  for  .".uggesting  such  a  “tripping”  device.  This  ring  device  should 
more  appropriately  be  named  “vortex  generator”  ring  due  to  its  major  effect  on  creating 
a  new  vortex-dominated  flow  field,  rather  than  to  its  effect  cn  “creating  separation”.  The 
usual  notion  of  “tripping”,  induced  from  the  laminar  to  turbulent  boundary  layer  tripping 
devices,  is  wrong.  However,  its  use,  unfortunately,  is  wide  spread  among  researchers. 

The  present  study  concerns  itself  with  projectiles  v;ith  “tripping”  rings  of  interest  to 
the  U.S.  Army.  In  particular,  the  interest  is  drawn  from  the  105  mm  M4S9  spike-nosed 
projectile.  This  round  is  the  training  round  for  the  105  mm  M456A1  fin-stabilized,  High 
Explosive  Anti-Tank  (HEAT)  projectile.  Computations  at  three  different  Mach  numbers  are 
performed  and  the  sting  support  effect  is  studied  to  assess  its  influence.  The  results  are 
compared  with  wind  tunnel  data  for  verification.  A  new  and  interesting  flow  pattern  was 
observed  and  is  described.  Also,  for  a  high  drag  case,  the  high  drag  solution  was  induced  by 
a  technique  described  earlier  in  Reference  5.  Further,  a  “buzzing”  flow  Ccise  was  encountered 
and  the  frequency  of  buzzing  Weis  computed. 
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II.  PROJECTILE  GEOMETRY  AND  THE  TEST  CASES 


The  projectile  configuration  of  interest  to  the  U.S.  Army  is  the  105  mm  M490  spin- 
stabilized  projectile  of  the  M489  training  round.  Wind  tunnel  test  results  were  available 
from  Reference  6  at  only  two  Mach  numbers;  3.0  and  3.5.  The  wind  tunnel  models  were  of 
2.0-inch  diameter  which  represent  only  48%  of  the  full  scale  projectile.  Several  variations  in 
the  projectile  spike  diameter,  projectile  shoulder  geometry,  and  the  afterbody  were  tested. 
The  chosen  configuration  for  the  present  study  was  model  number  1112,  which  is  depicted 
in  Figure  1.  The  tests  were  run  in  a  blow-down  supersonic  wind  tunnel  with  the  following 
test  conditions:  Mach  numbers  (3.0,3.5),  total  pressure  (60.0,70.0)  psi,  total  temperature 
(59.0,59.0)  and  Reynolds  number  (10.2,9.1)  xlO®  per  foot,  respectively.  Some  range 
firing  tests  were  also  made  for  the  spinning  full  scale  projectile^.  However,  due  to  the 
yawing  motion  observed  and  lack  of  documentation,  the  zero  yaw  axial  force  coefficient 
was  not  readily  available.  Also  it  should  be  noted  that  the  scaling  laws  between  the  full 
scale  and  the  0.48  scaled  models  for  spike-nosed  bodies  are  not  known.  The  flow  pattern 
and  the  flow  over  the  spike  are  highly  sensitive  to  the  spike  length  (in  association  with  the 
local  Reynolds  number  influencing  the  separation  region  along  the  spike)  and  in  conjunction 
with  the  projectile  shoulder  diameter  (the  shoulder  diam.eter  to  the  spike  diameter  ratio). 
Therefore,  the  wind  tunnel  results  were  favored  for  the  numerical  simulation. 

The  sting  support  was  0.67  inch  in  diameter,  corresponding  to  0.333  of  the  model  btise 
diameter.  Cases  of  angles  of  attack  between  ±8®  were  tested  and  measured.  No  error  bounds 
on  any  of  the  given  measurements  were  given  in  the  report  of  Reference  6. 


III.  COMPUTATIONAL  TECHNIQUE 

1.  GOVERNING  EQUATIONS 

The  compressible,  turbulent  Navier-Stokes  equations  for  axisymmetricand  two-dimensional 
flow  can  be  expressed®  in  the  following  strong  conservation  form  where  the  dependent  vari¬ 
ables  p,  u,  V,  and  e  are  mass  averaged  where  e  being  the  specific  total  energy,  u  and  v  are  the 
axial  and  radial  velocity  components,  p  and  p  being  mean  density  and  pressure,  respectively, 
and  t  is  time: 


where 


dQ'  dE'  OF' 
dt  ^  dx  ^  dy  ^ 


■  p 

■  pu 

pu 

puu  -^p-Txx 

pv 

• 

PUV  -  Txy 

.  . 

(pe  -i-  p)u  -  TxxU  -  TxyV  4-  qx  . 
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The  laminar  urbulent  Prandtl  numbers,  Pr  and  Prt,  were  assumed  constant  with 
values  of  0.72  anu  ,  respectively.  The  ratio  of  the  specific  heats,  7,  was  also  assumed 
constant  and  equal  to  1.4.  C\  and  Cp  are  specific  heat  caj^acities  at  constant  volume  and 
constant  pressure,  respectively. 

__  i*AAn 


=  4290 


sec*  — 


for  air.  The  total  energy  per  unit  mass,  c,  is  given  by: 

c  =  C„r  +  (l/2)(u"  +  v*). 

In  the  i  —  1)  computational  plane.  Equation  1  is  transformed  to  the  conservation  law 
form  and  the  equations  can  be  found,  for  example,  in  Reference  8. 


2.  TURBULENCE  MODEL 

Turbulence  is  modeled  through  a  modification  to  the  eddy  viscosity  model  of  Baldwin 
and  Lomax®.  This  widely  applied  model  employs  the  two-layer  concept  (inner  and  outer 
layers).  The  inner  layer  is  near  the  wall  and  is  modeled  as: 


€,•  = 


/  =  /:j/  1  —  exp 


(^)] 


The  magnitude  of  the  vorticity  |u;|  is: 


and  where 


I  I 


The  distance  normal  to  the  surface  is  y,  =  26,  k  =  0.40  is  the  von  Karman  constant,  and 
the  subscript  w  denotes  values  at  the  surface. 

The  model  switches  from  the  inner  to  the  outer  region  at  the  smallest  value  of  y  which 
the  inner  and  outer  values  cf  the  eddy  viscosity  are  equal  (i.e,  c,  =  Co).  The  e  for  the  outer 
layer  is  given  by: 

^0  ~  pP  CcpP naxy max PkLEB  (6) 
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A'  =  0.0168,  Cep  =  1.6,  Ckleb  -  0.3.  (9) 

Due  to  the  perpendicular  intersection  of  the  spike  surfaces  at  the  nose  tip  and  also  at  the 
facing  shoulder,  the  normal  distance  to  these  walls  “j/”  in  Equations  2-8  is  difficult  to  assign 
This  problem  was  approached  in  Reference  10  by  measuring  the  y  along  a  45“  ray  emanating 
from  the  point  of  intersection  of  the  two  perpendicular  walls. 


3.  THE  CODE 

The  computer  rode  was  developed  by  Patel  et  al®  and  utilizes  the  robust  explicit,  time 
dependent  method  of  McCormack,  The  code  is  vectorized  and  is  run  on  Cray-XMP/48 
machine.  However  the  present  computations  were  all  run  in  serial  arithmetic  mode.  The 
zonal  grid  and  overlapping  is  provided  in  the  code  as  represented  by  eight  different  zones  (can 
be  increased  if  so  desired)  The  user  prescribes  the  overlapping  between  regions  along  one 
line  of  adjacent  zones  (interface).  A  global  uniform  time  step  was  used  herein  rather  than 
grid-varying  time  steps  to  simulate  “time-accurate"  solutions.  The  time  step  is  determined 
from  the  CFL  (Courant-Fredrich-Levy)  condition,  with  a  factor  of  0.6  being  used  as  the 
Courant  number.  One  characteristic  time  length  for  the  flow  field  at  M  =3.0  and  3.5  was 
454.  xl0~*  and  389.  xlO"*  second,  respectively.  With  an  average  time  step  of  0.19  xl0“® 
sec,  one  characteristic  time  length  corresponds  to  2390  and  2050  numerical  time  steps  at 
Mach  =  3.0  and  3.5  ,  respectively  . 


4.  BOUNDARY  CONDITIONS 


No-slip  conditions  are  specified  on  all  wall  surfaces.  The  incoming  flow  conditions  are 
assumed  to  be  of  uniform  profiles  with  freestream  values  based  on  the  wind  tunnel  test 
values. 


A/oo 

roo(R“) 

Pooipsi) 

Rc{pcrft) 

1.9 

302 

7.31 

12.62  X  10® 

3.0 

185 

1.63 

10.20  X  10® 

3.5 

150 

0.98 

9.10  X  10® 

The  outgoing  flow  conditions  at  the  downstream  end  of  the  flow  field  were  imposed  as 
zero  gradients  along  the  body  axis  direction  at  a  distance  of  4.0  body  diameters  behind  the 
projectile  base. 
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The  outer  boundary  conditions  were  imposed  as  non  reflective  conditions,  i.e..  zero 
gradient  conditions  along  characteristic  lines  for  all  variables.  The  characteristic  direction 
is  determined  from  the  local  velocity  and  temperature.  This  approach  allows  setting  the 
“outer”  field  close  to  the  body  without  the  penalty  of  any  unnecessary  approximation  re¬ 
garding  shock  reflection,  or  degradation  due  to  far-field  conditions  (no-gradients)  being  set 
too  close  . 

At  the  symmetry  line,  ahead  of  the  spike  tip  and  behind  the  base,  a  two-point  zero 
gradient  boundary  condition  is  imposed  on  the  solved  variables. 


5.  INITIAL  CONDITIONS 

Computations  for  the  first  case  of  Mach  3.0  were  started  using  free  stream  values  ev¬ 
erywhere  in  the  domain  for  all  variables,  based  on  an  initial  lower  Mach  number  which  wais 
selected  as  1.5.  These  values  are  for  free  stream  pressure,  and  temperature.  The  density  and 
specific  total  energy  are  computed  accordingly,  using  the  equation  of  state  and  the  definition 
of  the  specific  total  energy.  The  rest  of  the  cases  were  started  in  the  same  manner  except  for 
the  case  of  Mach  3.5  with  no  sting,  which  was  started  using  the  solution  for  the  case  with 
the  sting  at  the  same  Mach  number. 

6.  THE  GRID 

Six  different  grid  zones  were  used  in  the  computation.  Those  zones  and  the  extent  of 
the  computational  domain  are  depicted  in  Figure  2. 

For  the  case  with  a  sting,  the  grid  used  for  the  six  zones  are  (20,85),  (13,72),  (7,65), 
(60,78),  (130,36),  and  (80,31)  respectively.  The  first  and  second  arguments  in  the  brackets 
refer  to  the  axial  and  radial  directions,  respectively.  This  grid  is  equivalent  to  16531  total 
points  or  about  a  (129x129)  grid.  For  the  case  of  no  sting,  the  same  grid  is  used  except  the 
last  zone  which  was  incre2ised  from  (80x31)  to  (80x51)  to  account  for  the  space  previously 
occupied  by  the  sting. 

One  restriction  in  the  present  grid  overlapping  technique  is  that  no  interpolation  is 
allowed  atl  the  interface  line  between  zones.  Thus,  each  point  on  either  side  of  any  two  zones 
must  have  the  exact  same  coordinates.  This  restriction  represents  some  constraint  in  the 
flexibility  if  the  grid  distribution  and  may  be  alleviated  in  future  development  of  the  code. 
Meanwhilej  to  accommodate  this  restriction  one  h£is  to  accept  unnecessary  clustering  of 
points  in  some  locations.  Figures  3  and  4  show  the  overall  grid  distribution  for  the  projectile 
configuratidn  with  and  without  sting,  respectively.  Figure  5  shows  the  details  of  the  grid 
near  the  spike,  depicting  the  clustered  points  along  lines  parallel  to  the  top  body  surface 
where  clustering  is  needed  near  the  body  to  resolve  the  turbulent  boundary  layer. 
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IV.  RESULTS 


1.  THE  CASE  AT  M  =  3.0 

This  case  was  started  using  the  earlier  stated  initial  conditions  and  the  previously 
described  grid  with  the  sting.  A  new  flow  field,  dominated  by  the  vortex  and  wake  generated 
by  the  “tripping”  ring,  extended  over  the  projectile  middle  body  and  beyond  the  bcise  of 
the  projectile.  This  vortex-wake  stre2Lm  interfered  with  the  usual  expansion  v'ave  expected 
at  the  base  for  the  supersonic  flow.  The  wake  flow  eliminated  the  expansion  wave  and  the 
“outer”  supersonic  flow  slipped  over  the  wake  stream  almost  unturned  thus,  no  expansion 
wave  took  place.  At  first,  this  behavior  was  suspect  but  was  later  confirmed  when  the  spark 
shadowgraph  photos  for  the  corresponding  full  scale  projectile  (obtained  during  the  range 
test  firing  of  Reference  7)  were  reviewed.  Examination  of  the  shadowgraphs  immediately 
confirmed  the  computations  and  showed  little  trace  of  the  expected  expansion  wave.  In 
addition,  the  shadowgraph  showed  very  clearly  the  vortex-wake  region  emanating  from  the 
vortex  ring  and  flowing  all  the  way  over  the  projectile  body.  This  flow  field,  now  appearing  to 
I  be  logical  and  acceptable  was,  nevertheless,  not  expected.  This  is  the  first  time  in  the  known 
literature  that  such  particular  flow  field  has  been  analyzed  and  computed.  A  schematic  of 
the  flow  pattern  is  shown  in  Figure  6  and  is  compared  with  the  typical  expected  flow  for  a 
spiked  body  without  a  tripping  ring.  This  vortex-dominated  field  effect  on  the  expansion 
f  wave  at  the  bzise  was  never  shown  in  the  literature  and  wcis  not  exposed  in  Mikhail’s  work^ 
since  no  base  region  computations  were  made  for  the  spiked  configurations  with  a  tripping 
[  ring. 

!  The  flow  field  is  depicted  in  Figures  7-9  by  the  Mach-line  contours  for  the  regions  near 

j  the  spike  tip,  the  shoulder,  and  the  base  of  the  projectile.  The  axial  force  coefficient,  here 
I  also  is  the  drag  coefficient,  was  computed  to  be  0.253  which  compares  extremely  well  with 
j  the  value  of  0.255  of  the  wind  tunnel  data®.  This  drag  comparison  is  depicted  in  Figure  10. 

The  case  took  about  40,000  time  steps  to  converge  only  because  there  is  a  very  small 
eddy  continuously  breaking  down  and  reforming  near  the  base  and  which  propagates  upwind 
to  the  frontal  spike  tip.  It  is  not  known  at  this  point  whether  this  unsteadiness  is  truely 
■  physical  or  purely  numerical  in  nature.  This  caused  the  drag  coefficient  to  vary  slightly  at  a 
very  slow  frequency.  The  drag  coefficient  variation  is  provided  in  Figure  11  as  a  drag  history 
plot. 

Figure  12  provides  the  anatomy  of  the  drag  for  the  case  at  M  =  3.0,  showing  the  value 
of  each  component  of  the  drag.  The  shoulder  drag  was  47%,  while  the  bcise  pressure  drag 
was  18%  of  the  total  drag.  The  spike  frontal  area  caused  30%  of  the  total  drag,  while  the 
skin  friction  of  the  whole  body  contributed  a  mere  0.2%.  The  vortex  ring  contributed  only 
4.4%. 

The  computation  time  on  the  Cray  XMP/48  was  about  40  minutes  for  each  1000  time 
steps.  For  a  40,000  steps  this  amounts  to  about  27  hours  of  CPU  time. 
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2.  THE  CASE  AT  M  =  3.5 


The  case  with  the  sting  was  started  from  the  same  uniform  initial  conditions  described 
earlier.  The  case  took  about  58,000  steps  to  converge  with  this  high  number  being  due,  as 
explained  above,  to  the  small  eddy  break-up.  The  computations  could  have  been  stopped 
at  heJf  the  number  of  steps  with  no  appreciable  difference.  However,  it  was  kept  running  to: 
first,  observe  any  changes;  and  second,  to  ensure  that  no  numerical  instability  would  occur. 

The  axial  force  coeflScient  computed  was  0.233  which,  as  expected,  is  lower  than  the 
value  for  the  M  =  3.0  case.  However,  the  tunnel  test  provided  a  higher  value  of  Co  =  0.309. 
It  was  immediately  realized  that  the  tunnel  test  resulted  in  the  high-drag  value  while  the 
computation  resulted  in  the,  low-drag  flow  pattern  value.  Also,  the  sting  was  suspected 
to  influence  the  frontal  spike  flow  pattern  through  this  oscillating  small  structure  eddy. 
Therefore,  the  effect  of  the  sting  had  first  to  be  assessed. 


a.  Effect  of  the  Sting 

The  sting  was  removed  and  the  new  grid  replacing  the  space  previously  occupied  by 
the  sting  was  constructed.  The  new  grid  is  depicted  in  Figure  4.  The  solution  was  started 
from  the  point  where  the  case  with  a  sting  ended.  The  computations  did  not  take  more 
thaii  10,000  steps  to  settle  again  and  converge.  Surprisingly,  no  appreciable  difference  was 
observed  in  the  frontal  spike  flow  or  in  the  flow  pattern.  Only  changes  in  the  base  flow  were 
observed,  as  expected.  Two  eddies  near  the  base  disappeared  when  the  sting  wais  removed. 
The  new  Cp  value  obtained  was  0.232  compared  to  the  0.233  for  the  case  with  the  sting. 
The  small  eddy  break-up  also  existed.  The  flow  field  for  the  case  with  and  without  sting  is 
given  in  Figures  13-14. 


b.  The  Higli  Drag  Case 

With  the  computed  axial  force  coefficient  being  lower  than  the  wind  tunnel  value, 
the  need  to  compute  the  corresponding  high  drag  ceise  was  necessary.  Beised  on  the  prior 
experience  cited  in  Reference  5,  the  high  drag  mode  was  then  “induced”  by  the  method  cited 
.which  suggested  a  pseudo  simulation  of  the  effect  of  an  angle  of  attack.  For  the  present  case, 
a  cross  flow  velocity  of  0.07  of  the  freestream  velocity  value  was  imposed  as  a  boundary 
condition  for  the  incoming  flow  (a  pseudo  4-degree  yaw  as  described  in  Reference  5).  The 
sting  was  put  back  for  the  simulation  of  the  existing  wind  tunnel  case.  The  computations 
were  started  from  the  previous  solution  of  M  =  3.5  case  with  a  sting. 

The  computation  took  30,000  steps  to  converge  (20  hours  CPU  time).  The  computed 
Cp  was  0.347.  This  value  obtained  falls  between  the  scatter  of  the  two  wind  tunnel  data 
points  of  0.309  and  0.3S6.  The  latter  value  was  for  an  identical  projectile  which  was  one  inch 
shorter  in  the  body  length  (near  the  base  of  the  projectile)  than  the  present  configuration. 

The  flow  fields  for  -both  the  low  and  high-drag  cases  are  very  similar  except  at  the 
base  near  the  sting.  The  flow  fields  are  depicted  in  Figures  15-16.  Figure  17  shows  the 
drag  coefficient  comparison  with  the  wind  tunnel  measurement.  The  numerical  convergence 
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history  is  given  in  Figure  18. 

3.  THE  CASE  AT  M  =  1.9 

After  the  computations  of  the  above  two  cases,  the  variation  of  Co  with  the  Marh 
number  was  of  interest  due  to  the  low-  and  high-  drag  cases  encountered  at  the  M  =  3.5 
case.  Although  no  wind  tunnel  data  is  available  at  lower  speeds,  an  additional  case  at 
A/  =  1.9  was  chosen  for  computation  to  assess  the  Cd  variation  with  Mach  number  for  this 
particular  projectile  configuration.  Having  no  wind  tunnel  data  to  compare  with,  and  after 
the  previous  result  that  the  sting  had  no  effect  on  the  spike  flow,  the  M  =  1.9  case  with  no 
sting  was  computed  to  simulate  the  in-flight  behavior  of  the  projectile. 

The  computation  started  from  the  same  uniform  initial  conditions  and,  after  a  large 
number  of  time  steps,  the  solution  showed  a  very  clear  type  of  oscillatory  (buzzing)  behavior. 
This  buzzing  is  of  a  different  class  than  the  small  unsteadiness  observed  earlier  due  to  the 
small  eddy  break-up.  The  clear  oscillatory  beh...'ior  was  observed  to  be  cyclic  after  about 
80,000  steps  and  persisted  to  110,000  steps  with  a  very  constant  frequency  of  about  3000 
steps  per  cycle.  Therefore,  ten  complete  cycles  were  encountered  before  the  computations 
were  stopped.  This  buzzing  flow  pattern  had  no  equal  amplitude  in  the  pressure  oscillation, 
attributable  to  the  non-cyclic  eddy  break-up.  This  same  phenomenon  is  frequently  observed 
in  combustor  flow  applications,  with  and  without  combustion  processes.  For  example  Jou  and 
Menon”,  in  their  study  of  non  reacting  flow  inside  a  ramjet  combustor,  obtained  a  buzzing 
flow  with  a  variable  pressure  amplitude  due  to  the  large  eddies  existing  in  the  flow  peist  the 
sudden  expansion  section.  The  flow  in  their  case  oscillated  cyclically  with  a  favored  frequency 
but  with  variable  pressure  amplitude.  The  same  phenomenon  there  was  also  explained  as 
the  effect  of  continuous  formation  and  break-up  of  large  eddies  in  the  flow. 

For  the  present  case,  with  the  time  step  At  =  0.19  x  10“®  sec,  the  frequency  of  the  buzz 
was  5333  Hz.  The  drag  coefficient  oscillated  between  about  0.22  and  0.62  with  a  “mean” 
value  of  about  0.42.  The  flow  pulsated  between  high  pressure  and  low  pressure  points. 
The  oscillatory  drag  behavior  is  depicted  in  Figure  19  and  the  more  detailed  oscillatory 
cycle  is  shown  in  Figure  20.  The  unequal  amplitude  of  the  oscillation  was  disturbing  at 
first,  but  when  consulting  the  literature,  the  same  behavior  was  readily  confirmed.  Jou 
and  Menon”obtained  a  large  eddy  break-up  pattern  when  computing  non-reacting  flow  in  a 
dump  combustor.  The  flow  in  their  case  also  oscillated  with  non-equal  amplitude. 

The  flow  field  for  this  interesting  buzzing  flow  is  captured  both  at  the  high  and  low 
pressure  points  of  the  cycle.  Figures  21-26  depict  the  flow  near  the  spike  tip,  the  shoulder, 
and  the  projectile  base.  Some  flow  changes  are  observed  as  the  flow  pulsates  very  rapidly. 

Finally,  the  drag  coefficient  variation  with  the  Mach  number  is  depicted  in  Figure  27  for 
all  cases  computed  and  is  compared  with  the  available  wind  tunnel  data.  The  data  validates 
the  computations  at  Mach  3.0  and  3.5,  while  the  computations  extend  the  results  down  to 
Mach  1.9. 
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V.  SUMMARY  AND  CONCLUSIONS 


Five  computations  were  made  for  a  spike-nosed  projectile  with  a  vortex  ring  at  Mach  1 .9, 
3.0,  eind  3.5  at  zero  angle  of  attack.  The  computations  were  made  by  solving  Navier-Stokes 
equations  in  an  overlapped  zonal  gridding  topology,  using  McCormack’s  explicit  scheme. 
The  results  obtained  are  summarized  as  follows: 

1.  A  new  flow  pattern  which  is  dominated  by  the  vortex  and  wake  generated  by  the 
“tripping”  ring  is  computed  and  described  .  The  vortex-wake  field  extends  to  the  base  of 
the  projectile  and  interfers  with  the  usual  flow  expansion  there.  This  flow  was  verified  by 
very  clecir  firing  range  spark  shadowgraph  photos  for  the  full  scale  model  tests^.  The  axial 
force  coeffidient  w^  in  excellent  agreement  with  the  wind  tunnel  data  at  Mach  3.0. 

2.  At  Mach  3.5,  both  low  and  high  drag  patterns  were  obtained.  The  high  drag  value 
agrees  with  the  wind  tunnel  data.  This  experimental  data  was  reported'  without  explanation 
of  why  the  drag  was  higher  at  Mach  3.5  than  at  Mach  3.0. 

3.  The  sting  effect  on  the  flow  pattern  and  axial  force  wais  studied  and  found  to  be 
almost  non-existent  for  the  present  sting  diameter  which  was  only  0.33  of  the  projectile 
dizuneter.  The  base  flow  pattern,  however,  was  influenced  and  showed  at  least  two  eddies, 
which  vanished  when  the  sting  was  removed,  for  the  case  of  Mach  3.5. 

4.  A  detailed  drag  anatomy  was  provided  ,  shedding  a  very  interesting  light  onto  the 
total  drag  picture  of  the  projectile.  The  front  spike  contributed  30%,  while  the  “shoulder” 
drag  is  47%  of  the  total  drag  for  the  case  at  Mach  3.0. 

5.  At  the  lower  Mach  number  of  1.9,  the  interesting  “buzzing”  (unsteady)  flow  wcis 
encountered  and  was  computed  successfully.  The  frequency  of  oscillation  was  found  to  be 
5333  Hz.  The  axial  force  coefficient  oscillated  between  0.22  and  0.62,  with  a  “mean”  value  of 
0.42.  Mo  wind  tunnel  tests  were  done  at  that  Mach  number.  Therefore,  the  designers  of  the 
projectile  were  apparently  not  aware  of  the  possible  buzzing  flow  at  lower  Mach  numbers. 

The  numerical  procedure  proved  to  be  robust  and  reliable  once  “enough”  grid  points 
are  used.  The  computations  proceeded  in  a  straight  forward  manner  and  without  major 
interruptions,  adding  to  the  previous  success  reported  in  Reference  5  for  spike-nosed  bodies. 
Possible  future  improvements  should  consider  better  turbulence  model  for  the  base  region, 
as  well  as  more  thorough  modeling  for  turbulence  in  the  recirculating  flow  region  ahead  of 
the  projectile  shoulder. 

Finally,  it  is  recommended  to  further  test  the  present  numerical  procedure  on  more 
complex  spike-nosed  configurations,  such  as  for  those  with  fin  boom,  boattail,  and  fins. 
Eflforts  are  currently  being  made  to  follow  up  on  this  recommendation. 
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e  5.  Details  of  the  grid  near  the  spike 
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Figure  8.  Flow  field  for  M  =  3.0  case  -  Mach-line  contours  near  projectile  shoulder 
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Figure  9.  Flow  field  for  .M  =  3.0  case  -  Mach-line  contours  near  the  sting  support 
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Figure  10. Drag  coefficient  comparison  for  M  =  3.0  case 
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Figure  14.  Flow  field  for  M  =  3.5  case  -  Mach-line  contours  for  case  without  sting 


Figure  20.  Oscillatory  cycle  for  the  M  —  1.9  case 


Figure  25.  Flow  field  for  the  M  =  1.9  case  -  near  the  base,  at  high  pressure  point 
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LIST  OF  SYMBOLS 


Art}  =  reference  area,  (7rcP/4) 

Cd  =  drag  coefBcient,  drag  force/(0.5/)V^i4re/) 

Cp  =  specific  heat  under  constant  pressure 

Cv  =  specific  heat  under  constant  volume 

d  =  reference  diameter 

ds  =  spike  diameter 

e  =  specific  total  energy 

J  =  Jacobian  of  the  coordinate  transformation 

M  =  Mach  number 

p  =  static  pressure 

Re  —  Reynolds  number 

T  =  temperature 

u,v  =  velocity  components  in  the  X,  y  directions 

Voo  =  free  stream  velocity 

X,  y  =r  Cartesian  coordinates  for  2-D  ceise,  axial  and  radial  coordinates  for 

axisymmetric  case 


Greek  Symbols 
o  =  angle  of  attack 

7  =  ratio  of  specific  heats 

p  =  density 

p  =  laminar  (molecular)  viscosity  coefficient 

e  =  turbulent  eddy  viscosity  coefficient 

(tV  ~  transformed  coordinates  in  the  computational  plane  for  the  coordinates  x,y 


00 


=  denotes  total  (stagnation)  condition 
=  free  stream  condition 
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